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I. INTRODUCTION 


A. Project Objectives and Approach 

The overall objective of the National Park Service (NPS) 
Applications System Verification Test (ASVT) Project is to develop 
the capability within the NPS to use computer-implemented techni- 
ques to derive geographically referenced surface cover information 
from Landsat Multispectral Scanner (MSS) data. For purposes 
of this project, techniques were to be developed and tested in 
three national park areas: 

- Olympic National Park 

- Death Valley National Monument 

- Shenandoah National Park 

Phase I of the NPS ASVT project required the processing 
and classification of four seasonal Landsat data sets for each 
park using National Aeronautics and Space Administration (NASA) 
developed computer programs; selection of one data set or the 
mergence of two or more data sets to improve surface cover class- 
ifications; data base building (where appropriate) for input of 
topographic data, etc. as necessary; accuracy verification; 
map product preparation; NPS evaluation; and service-wide review 
for technology transfer. 

Phase II is to be initiated after successful evaluation 
of products for each of the three parks by NPS personnel. This 
phase involves NPS acquisition of equipment (computer hardware 
and software), NPS staff recruitment and hiring, training of NPS 
personnel in use of NASA programs by NASA personnel at the National 
Space Technology Laboratories (NSTL) , Earth Resources Laboratory (ERL) , 
and start-up of the .system under ERL guidance. 
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This report details the results of Phase I of the 
cooperative NPS/ASVT project with respect to Shenandoah National 
Park (see also, Cibula 1980 and 1981). 

B. Description of Shenandoah National Park 

Shenandoah National Park lies within the northern portion 
of the Blue Ridge Mountains, 70 miles west of Washington, DC (Fig- 
ure 1) . The Shenandoah River flows through the Shenandoah valley 
which is located to the west of the park. This valley is flanked 
by the Blue Ridge Mountains on the east and the Allegheny Moun- 
tians to the west. The forty mile long Massanutten Mountains di- 
vide part of this valley and separate the north and south forks of 
the Shenandoah River. The Piedmont Plateau lies to the east of 
the park. 

From the mid-eighteenth century and into the early 1900 's 
the hollows, slopes and ridges of the lands which now constitute 
the park, were the home of mountain people. These pioneers existed 
marginally by growing vegetables and apples, raising hogs, chickens 
and cows, utilizing nuts from hardwoods in the forest (which then 
included the American Chestnut) , as well as harvesting wild honey 
and trapping for furs. The production of moonshine played an im- 
portant role in their economy. 

Early in this century, these people began to abandon this 
land as the soil steadily became less productive. Lumber companies 
had depleted the forests, most game was gone, the soil was exhausted 
and, as the result of misuse, much of it had eroded down the slopes. 
Thu American Chestnut ( Castanea dentata ) was one of the dominant 
species of these forests and in most years provided a bountiful 
crop of nuts. During the late 1920's, the chestnut blight ( Endothia 
parasitica ) , eliminated this tree from these forests, thereby 
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completely destroying a staple of the local economy. 

As a result of this earlier exploitation, the forest 
cover is essentially second growth mixed hardwoods with various 
oaks predominating. Large areas of meadowland existed earlier 
as a result of clearing for croplands and grazing. Most of 
these have reverted to forest, being covered now with various 
successional tree and shrub species. 

Red oak ( Quercus rubra ) , dominates the higher elevation 
forests and is one of the last hardwoods to leaf out in spring. 

At high elevations, stands of this oak may be nearly pure, but 
at lowe.' elevations it is found mixed with Chestnut oak (^. prinus ) , 
white oak (^. alba ) and the hickories ( Carya spp . ) . Near the 
summit of Hawks Bill Mountain, red oak is associated with red 
spruce ( Picea rubens ) and balsam fir ( Abies balsamea ) where these 
latter two species persi.^t as relics from the Pleistocene Period. 
Within the park, the red oak forest is the most prevalent, occuring 
over a fairly broad range of environmental conditions since this 
forest occupies the more mesic ridge tops as well as the upper side 
slopes. Often these forests have an understory of mountain laurel 
( Kalmia latifolia ) which is most commonly found on dry rocky areas 
and ridges beneath the oak. 

The second most prevalent forest type is the chestnut oak. 
This forest occupies the slightly drier side slopes with a southern 
or southwest aspect, and is dominated by chestnut oak with red oak 
as its primary associate. With the removal of the American Chestnut, 
the expected climax forest under present conditions would probably 
be oak/hickory (Braun, 1950). 
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White oak/Scarlet oak (^. cocctnea ) aie more predominant 
at lower elevations and often are associated with chestnut oak, 
pignut hickory ( Carya glabra ) , scrub or bear oak (^. ilicifolia ) . 

Cove hardwoods comprise a mixture of moist site species, 
and includes black birch ( Betula lenta ) , yellow or tulip poplar 
( Liriodendron tulipifera ) , red maple ( Acer rubrum ) , various eims 
( Ulmus spp.), and pawpaw ( Asimina triloba ) . In some areas, tulip 
poplar forms dominant stands where this species represents one 
of the pioneer species in reclaiming to forest old abandoned 
fields or pastures. The cove hardwood forest is the richest and 
most diverse forest type. 

Various species of pine occur in several different ecolo- 
gical regimes within the park. Hemlock ( Tsuga canadensis ) stands 
are few and scattered, the best known being the Limberlost area 
which represents ore of the few areas of virgin forest within 
Shenandoah National Park. Hemlock is most commonly found along 
streams at their upper reaches at higher elevations within the 
park. Common associates with hemlock are red maple, black and 
yellow birch, red oak and white ash ( Fraxinus americana ) . 

White pine ( Pinus '’ trobus ) is most commonly found on more 
xeric sites than hemlock, and usually is fo’md in association with 
red, chestnut and white oaks, flowering dogwood ( Coi as florida ) , 
black locust ( Robinia pseudoacacia ) , white ash, and sassafras 
( Sassafras albidum ) . 

Virginia pine (Pi nus virginiana ) , occurs on former fields 
and dry sites below 2500' elevation. This pine is early success- 
ional but does not persist as long as white pine due to its shorter 
life. 
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All pines are early successlonal as they are 
Intolerant of shade and must germinate and establish seedlings 
in relatively open situations such as abandoned fields or areas 
in the forest where the canopy has been opened by some sort of 
catastrophic occurance. 

From the above, and the information available on the 
earlier history of these lands, it is evident that most of the 
forest cf Shenandoah National Park is relativ ly new, as a re- 
placement for many abandoned fields and farms. As such, this 
forest is in an early stage of succession , and cannot be expected 
to achieve climax conditions until sometime in the relatively 
distant future. The successlonal aspect of these forests make 
them more difficult to describe and typify. 

II. METHODOLOGY 

A. Data Acquisition 

1 Landsat frames chosen for study: 

Landsat satellites have been collecting data over the 
earth's surface since 1972. The satellities are in near polar, 
circula::, sun-synchronous orbits with an altitude above the 
earth's surface of approximately 920 km (570 miles), and have 
a nominal 9:30 am crossing of the equator during the descending 
node. They circle the earth every 103 minutes (14 times a day) 
with each successive pass displaced to the west approximately 
26° of longitude due to the earth's rotation. The multispectral 
scanner (MSS) , the primary sensor aboard Landsat provides a 
continuous image of a strip of the earth's surface 185 km (115 
statute miles) wide. On the fifteenth pass, occurring 24 hours 
after the first pass, the coverage is shifted to the west an 
amount that provides a sidelap with the previous day coverage of 
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147o at the equator to 1007o at the poles. After 18 days, orbit 
252 (on the 19th day) retraces that of the first orbit providing 
repetitive coverage (U.S. Geological Survey, 1979). 

Data analyzed in this study were obtained from the MSS. 

The MSS measures radiance in the following wave length bands : 
0.5-0.6uni (green), 0.6-0. 7ym (red), 0.7-0.8ym and 0.8-1. lym 
(both near infrared) . Values obtained in each of these bands 
for each ground resolution element (pixel) form a multispectral 
data set which was the basis for analysis. The four radiance 
values for each pixel form a "data vector" which is associated 
with that element. Use of a digital computer permits each data 
vector to be assigned to a class which contain vectors of similar 
character, hence a classification. This aspect will be discusses 
in more detail in the appendix which details the data processing. 

The analog signals produced by each of the four MSS detectors 
are digitized and formatted into a 15 megabit data stream for 
transmission to an Earth receiving station. After being received 
by the NASA Image Processing Facility (IPF) , the MSS data are 
transformed into segmented imagery with an overlap between "frames" 
(frame size 185 km x 185 km) . The framed data are also produced 
onto 9 track 800 or 1600-bpi computer compatible tapes (CCT's). 
These are the tapes which are available on order from the USGS 
EROS Data Center, Sioux Falls, South Dakota and which were used 
in this study. 
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Four Landsat frames, corresponding to one of the 
four seasons, were chosen for the Shenandoah National Park study. 
These were : 

- 9 May 1976; 2473-15081; Chosen to provide cover'Jge 
during spring when most broadleafed hardwoods would be expected 
to leaf out. 

- 5 June 1976; 5413-14404: Selected to provide summer 

coverage, with the time chosen to be early in the season so 
most broadleafed trees would have mature but no physiologically 
aged leaves. 

- 18 October 1976; 2635-15034: Selected so that at this 

date, peak fall coloration would bt' achieved over most areas with- 
in the park. 

- 27 February 1976; 2401-15102: Chosen to provide coverage 

of a winter but snowless condition within the park. The date 
choice v;as made to give a higher sun angle than would have been 
available in December, thereby minimizing problems with shadows. 

2. Field studies 

After the Landsat frames described above were spec- 
trally classified (discussed in RESULTS section), selected areas 
representing each spectral class were chosen from the Landsat data 
and transferred to USGS 1:62,500 topographic maps for field use. 

In gathering and recording the field data, a modified 
Braun-Blanquet (1^51) technique was employed. The technique 
employed is discussed in more detail by Cibula and Nyquist (1981). 
Most sites were visited by ground field teams. A few of the 
more inaccessable sites were visited by a helicopter team. The 
canopy analysis employed permitted the stratification of forest 
sites, ususally into four distinct canopy levels. This is shown 
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schematically in Figure 2. An example of a completed ground 
truth form following these criteria, is shown as Figure 3. 

After completion of this form, a photograph 
was taken of each site using Kodacolor film. Later, after 
processing, the print (s) which corresponded to each site were 
mounted with the respective ground data sheet. The original 
data sheets are archived at the National Park Service Denver 
Service Center and Headquarters, Shenandoah National Park. 

Appendix III gives the UTM coordinates which define the center 
point for each site. 

This ground truth data became the basis by 
which the spectral classes were identified. Furthermore, as it 
evolved later, this data was invaluable in developing the 
relationships used to refine this classification through the 
use of topographic digital data. These concepts will be detailed 
later in this report. 

B. Data Processing 

Final processing and manipulation was accomplished with 
the Earth Resources Laboratory Applications Software (ELAS) . Al- 
though initial processing was undertaken before ELAS was completed, 
these earlier programs are now incorporated as ELAS overlays 
written as FORTRAN modules. As a result of the above and the 
fact that ELAS represents the software package that will be trans- 
ferred to a potential user of Landsat data (e.g. , the National 
Park Service) , the processing flow of Landsat data for Shenandoah 
National Park will be described in terms of ELAS overlays. 
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ELAS is divided into two major components, (a) the oper- 
ating subsystem and (b) the application modules. The operating 
subsystem is written in FORTRAN as are the application modules; 
application modules rely on the operating subsystem for machine- 
dependent functions. The EL/iS subsystem requires an interactive 
display system. In use with this project, ELAS operated on an 
Interdata (Perkin-Elmer) 8/32 computer with a Comtal series 8000 
image display system. The various overlays used in this project 
are described briefly in Appendix I. 

III. RESULTS 

A. Discussion of Classifications and Rationale for Final 

Choice 

All data after spectral classification by SRCH were 
screened by both NFS and NASA personnel. The spring data set 
(Frame 2473-15081) presented a distinct zonate pattern with a 
number of separate classes that appeared to be related to ele- 
vation. This aspect was not as evident on the other scenes 
classified. Further analysis revealed that this pattern was 
related to the different degrees of leaf emergence found at 
different elevations within the park. Because of the phono- 
logical aspect of this data set, it was decided not to consider 
this data further. The fall data set exhibited much the same 
pattern but to a lesser degree. In this case, it was the fall 
coloration which exhibited differences due to elevation. It 
was also decided not to use this data. 

The summer data set, Frame 5413-14404, did not exhibit 
these anomalies, rather, the differences in spectral classes were 
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either related to true differences in forest type, or as a 
result of differing slopes and aspects of the forest communi- 
ties. The major deficiency appeared to be inadequate spectral 
signatures for the various confier communities. This deficiency 
was overcome with the winter data (Frame 2401-15102) as there 
were 5 spectral signatures that appeared to be related to conifer 
or evergreen associations. In addition, it was learned that 
winter data could be used to detect the presence or absence of 
evergreen species in the lower canopies of the various hardwood 
communities since during winter, understory foliage was not ob- 
scured by the leafless upper canopy. With these observations, it 
was decided to use both the summer and winter scenes, and to merge 
these classes in the ruanner described below to obtain the final 
classification for Shenandoah National Park. 

Initial analysis was accomplished with the summer data se':. 
Figure 4 shows the two-space plot for this data. This plot is 
the average of the two visible channels (1 & 2) as the abscissa 
and the two infrared channels (3 & 4) as the ordinate. Table I 
gives the means for each class. 

Analysis of field data permitted a description of the land 
cover types to be written for each spectral class. During this 
analysis, the data base channel with the elevational data was also 
utilized.* This allowed a discrimination to be made in a number 
of instances where the land cover type at one elevation interval 
was different from the type found at another interval. Ir. practice, 
use of the elevational data permitted a specific elevational contour 


*See Appendix I 


pg. 40, and discussion below. 



TABLE I 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Means and Color Values for SRCH Classes on Frame 5413-14404 
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to be displayed in the classified data. These determinations 
were made with NFS representatives from both Shenandoah National 
Park and the Denver Service Center while viewing tne displayed 
data at NSTL. The result of these studies are given in Table II. 

This data suggests that a number of spectral classes which 
are heterogenous with respect to land cover could be resolved if 
an elevational division were involved. To accomplish this, the 
Programmable Calculator (PCAL) of ELAS was used (for a discussion 
of this and other ELAS options, see Appendix I). 

As an example, analysis of field data and the spectral plots 
(Figure 4) indicated that the red oak communities even at the 
date of this Landsat pass (June 5), had not fully leafed out. 
Nearly all other deciduous forest species were in full leaf this 
late in the season. It is known that among hardwoods, red oak is 
about the last to leaf out (Fowells, 1965). This phenology allows 
us to separate this forest community from others which are already 
in full leaf. 

Continuing with this example, examination of Table II demon- 
strates that spectral classes 11, 15, 16, 17, 18, & 19 represent 
red oak at elevations above 2500'. These same spectral classes 
represent pasture or inert areas below this elevation. By use of 
the data base, elevational data is brought into correspondence 
with the classification. This means that for each and every 50m 
classification pixel, there exists a corresponding 50ir pixel in 
the elevational data where the numeric value of the elevation is 
the average for the terrain within this pixel. 
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TABLE II 


ANALYSIS OF SPECTRAL CLASSES 
Frame 5413-14404, from SRCH 


Spectral 

Class Description 

1 Grasses and clover; both low and high elevations, 
separate at 2,300 ft elevation 

2 Red oak dominant, SE to S aspect, slope 10° to 22° 

3 Primarily red and chestnut oak; some sites may have 
locust interspersed 

4 Oak/yellow poplar /hickory , NW & W aspect 

5 Red oak dominant with hickory; S and Fi-J aspect 

6 Chestnut oak/red oak; NW aspect 

7 Transitional between red oak communities which have 

not fully leafed out and lower elevational oak com- 
munities which are fully leafed; primarily composed 
of chestnut oak/red oak 

8 Grasses, unimproved pr.sture, old field succession at 
low elevation; high elevation meadows within park, 
separate at 2500 ft. 

9 Open canopy forest with grass understory; includes 
apple orchards 

10 Chestnut oak/red oak, S , W & NW aspect 

11 Low elevation pastures; at higher elevations is green 
underctory in oaks which have not leafed out. Separate 
at 2500 ft. 

12 Transitional between non-leafed oak and leafed oak. 
Above 1600 ft, this class represents red oak, while 
below this elevation, this class is representative 
of the red oak/chestnut oak association 

13 Oak/hickory /yellow poplar; E to SE aspect, slope 
to 50° 

14 Chestnut oak/red oak 

15 Mixed hardwoods; above 2500 ft, this represents red 
oak, below is representative of partially bare (inert 
mtls) areas. Separate at 2500 ft. 
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16 Similar to (15); separate at 2500 ft. 

17 Red oak above 2500 ft; pastures in low elevations, 
separate at 2500 ft. 

18 Red oak above 2500 ft; pasture at 1000 ft and lower; 
separate at 2500 ft. 

19 Pasture 500 ft to 1000 ft; above 2500 ft represents red 
oak. Separate at 2500 ft 

20 Red oak above 2500 ft 

21 Chestnut oak/red oak; W-NW aspect, 20° to 40° slope 

22 Red oak forests 

23 High reflecting inert materials 

24 Red oak forests at high elevations 

25 Red oak forests at high elevations 

26 Red oak forests at high elevations 

27 Red oak forests at high elevations 

28 Below 1000 ft lowland meadows, orchards or pasture lands; 
abo\e 2000 ft upland meadows 

29 Same as (28) above 

30 Dark inert, cultivated fields, roadbeds; etc., below 2500 
ft.; red oak forests above 2500 ft 

31 Dark inert below 2500 ft.; upland meadows above 2500 ft 

32 Mixed hardwoods, primarily red oak 

33 Very limited data on this class; appears to be regeneration 
over djsturbgd areas, correlates highly with east aspect 
and 25° - 35 slope 

34 Inert materials of very low reflectance 

35 Inert materials of moderately low reflectance 
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The use of che Programmable Calculator (PCAL) option in ELAS 
would permit a given class (of the original 35) to be split into 
two (or more) new clar.ses with the division based on elevation. 

In this example, class 11 remafned class 11 below 2500'. Above 
this elevation, this spectral class was assigned to class 37 
(remember, that the original SRCH class ' fication consi:»ted of 
35 spectral classes, while the computer system used has a capa- 
bility of handling 63 classes). Now, class 11 could be accu- 
rately described as lowland pastures while class 32 could 
accurately be described as part of the red oak forest community. 

In a similar manner, above 2500', classes 15, 16, 17, 18, and 19 
became 39, 40, 41, 42, and 43. These latter, now also could be 
described as components of the red oak community. 

In addition, PCAL could be used to combine specific classes 
of one classification with another classification to modify and 
expand the detail in the final product. In the winter data, 
pine was accurately pinpointed. In the summer data, tb" separation 
was not nearly as clear. As a result, it was desirable to write 
in over the summer data the known pine classes. In all, there 
were 13 separate PCAL operations performed which yielded the final 
classification shown, in Table III, Figures 5A and 5B. Details 
on each of the PCAL operations employed are given in Appendix II . 

B . A creage 

With the park boundary imbedded in the data, it was pos- 
sible to compute acreages for each of the land cover types which 
occurred within this boundary polygon. This gives directly, the 
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TABLE III 


ORIGINAL PAGE 13 
OF POOR QUALITY 


SHENANDOAH NATIONAL PARK 
Land Cover Classification 


Group 

Classes Included 

Description 

1 

1.57 

Grasses 6c forbes below level 2250 ft. 

2 

2,3,5,15,20.22,24, 

25,26.27.32,37,38, 

39,40,41,42,43,44, 

45.46,47.48,49,50 

Red Oak forests 

3 

4,13 

Yellow poplar; oak/hickory 

4 

6,7,9,10,12,14,21 

Chestnut oak/red oak 

5 

17,19,56 

Moderately bright inert materials 

6 

8 

Inert materials of low reflectance 

7 

11,30,31 

Inert materials of moderately low 
reflectance 

8 

18 

Pine 

9 

23 

Inert materials of high reflectance 

10 

33 

Vigorous forest regeneration, includes 
yellow poplar and other hardwoods that 
are not oak 

11 

34 

Inert materials of very low reflectanc 

12 

36,51,53 + (28. 

29 6c 32 above 2000') 

Upland meadows (Big Meadows) 

13 

52 

Defoliation in deciduous forests 

14 

54,60 

Hardwood forests (primarily red oak) 
with evergreen understory of Kalmia 
latifolia and Pinus strobus with some 
mixed grasses in some areas 

15 

55 

Red oak with evergreen understory 

16 

58,63 

Dense pine 

17 

59 

Mixed pine /hardwood 

18 

61 

Pine 
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Group 

Classes Included 

Description 


19 

62 

Pine alone or pine with a 
mixture of hardwoods 

slight 

20 

28 & 29 below 
2000' elevation 

Lowland meadow or pasture 

lands 

21 

Overlay of winter 
classes 4,6,8,10, 
20 above 2000' 

Upland meadows which have 
scattered woody perennials 
of the meadow community 

some small 
as part 

22 

8 

Old field succession 
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SHENANDOAH NATIONAL PAIiK 

Land Cover Classification 

•. * 

Descript ion 

Grasses & forbs below 2250’ 

Red Oak forests 

Yellow poplaf; oak/hickory 
» 

chestnut oak/red oak , , 

Moderately bright inert 

* materials 

Inert materials of low*^ 
reflectance i 

■ « 

Inert materiali^ of moderatley 
low reflectance • 

Pine 

Inert material's of high 
reflectance 

. Vigorous forest regeneration, 

. includes yellow poplar and other 
hardwoods that are not oak. 

Inert materials of very low » 
reflectance 

Upland meadows (Big Meadows) 

Defoliation in deciduous 
forests 

Hardwood forests (primarily red 
oak) with evjergreen ’ understory ^ 
of Kalmia latifolia and Pinus 
srro bus wit^i some mixed grasses 
. in some areas , * • • 

Red oak with evergreen under- 
, story • 

♦ «Dertse pine 

• Mixed pine /hardwood 
Pine 

Pine alone or pine with a 
sLi«ht mixture of hardwoods 

Lowland meadow or pasture 
lands 

Upland meadows yhich haA^e some ® 
small, scattered woody perennials 
as part of the. meadow community 

Old field succe&sion 


- it's 
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acreages of each land cover type within the park. In addition 
to the central area of the park, there existed also, four small 
separate islands. Each was calculated separately. These acreages 
are shown in Table IV. Group 60 is the embedded boundary. 

C . Accuracy 

For field verification, 200 northings and eastings were 
randomly generated for areas within the park. For each site, a 
15x15 matrix was printed with numbered groupings corresponding to 
the land cover groupings given in Table III. Each site was also 
plotted onto 1:62.500 USGS maps. The UTM coordinates for these 
sites are given in Appendix IV. 

In October 1079, field teams visited most of these sites 
and completed a canopy analysis as descirbed earlier, for the 
field sites relating to the spectral classes. Where possible, 
a photograph was taken of each site. Most sites were visited 
by helicopter, with a few of the more accessible sites being 
visited by ground teams. After collation of this data, each 
site was examined with respect to the land cover groupings given 
in Table III. A judgement was made from the field data as to 
which descriptive term best applied to the site. If the site was 
located at the junction of different land cover types, each was 
examined and assigned a group description. This completed, the 
matrices were examined in a 3x3 array about the center point (150 
meter square area) . The group number for the center point was 
noted, as well as any other group numbers which occurred within 
the 9x9 arrav. If there was coincidence between the group numbers 
obtained on the matrix printout and the group description ascribed 
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TABLE IV 


CLASSIFICATION ACREAGES FOR SHENANDOAH 
NATIONAL PARK 


Group 

Central 

Area 

1 

2 

3 

4 

Total 

1 

29 

5 

- 

- 

- 

34 

2 

81,464 

35 

106 

2 

9 

81,616 

3 

25,894 

7 

12 

4 

46 

25,963 

4 

32,936 

22 

2 

14 

85 

33,059 

5 

205 

20 

- 

- 

- 

225 

6 

162 

11 

- 

- 

- 

173 

7 

3,658 

1 

- 

- 

5 

3,664 

8 

156 

4 

- 

- 

- 

160 

9 

11 

22 

- 

- 

- 

33 

10 

1,636 

2 

1 

- 

- 

1,639 

11 

1,859 

- 

- 

- 

- 

1,859 

12 

566 

- 

- 

- 

- 

566 

13 

153 

- 

- 

- 

- 

153 

14 

13,130 

7 

33 

4 

3 

13,177 

15 

1,433 

17 

6 

- 

- 

1,455 

16 

3,244 

26 

2 

- 

2 

3,274 

17 

740 

16 

- 

- 

1 

756 

18 

1,183 

35 

- 

- 

1 

1,218 

19 

2,896 

7 

6 

- 

1 

2,910 

20 

83 

31 

- 

- 

- 

114 

21 

105 

- 

- 

- 

- 

105 

22 

43 

9 

- 

- 

- 

52 

60 

3,827 

36 

27 

13 

37 

3,940 

Total 

175,413 

311 

194 

38 

189 

176,145 


Total acres , 
whole park 

O *7 



to the site on the basis of field analysis, the site was considered 
correct. If correspondence did not exist, the group classification 
was considered to be in error. These results are summarized in 
Table V. 

Since time constraints dictated the use of helicopter teams, 
it was not always possible to determine if some forest communities 
had an evergreen understory as the trees were foliated at the time 
of the accuracy assessment. In these cases, where the classifi- 
cation group indicated red oak forests with an understory of 
Kalmia latifolia , this group was considered correct if the field 
data noted that red oak was present. 

Some field data, indicated a forest type not separated in the 
classification - black locust woodlands ( Robinia pseudoacacia ) . 

In most cases, these areas were classified as red oak communities. 

In these instances, the classification was considered incorrect, 
but it is significant to note that in almost every field site 
observed, red oak to some degree along with black locust comprised 
the community. 

Even with these variances, the overall accuracy of this clas- 
sification was 87.97o. 

IV. CONCLUSION 

This classification represents one of the first attempts to 
apply elevational data to several different seasonal Landsat frames 
with the final classification being the result of combining these 
completed classifications using a rigorously defined logic developed 
from the field work. In this classification, several features of 
use in park management resulted from this process. High elevation 
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meadows (both with and without recently emergent perennial 
vegetation) were separated, as well as areas in oak forests which 
have an evergreen understory as opposed to other areas which do 
not. This latter feature is of importance in wildlife habitat 
analysis . 
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APPENDIX I 


The various overlays available in ELAS and used in this 
project are descrioed briefly below. For more detailed infor- 
mation about these and other overlays and the ELAS software, 
the reader is referred to NASA/ERL (1980). 

REFO - Reformat of Goddard Landsat CCTs 

This module is used to reformat data in the Goddard tape 
format (Goddard Space Flight Center, 1976) into the ELAS data 
file format. This program however, does not reformat EROS com- 
puter compatible tapes (CCTs). Program NCCT (below) should be 
used for this purpose. 

NCCT - Reformat of EROS' CCTs 

This module reformats the EROS CCT tapes (EROS Data Center, 
1978) and writes the data in the standard ELAS data file format. 

JCOP - Join and Copy Data Files 

This module allows one or more files to be input and joined 
according to line and clement numbers. As an example, a Landsat 
scene is normally contained in four different files. Here JCOP 
can be used to put the entire scene or a subarea of the scene 
together in one fiie. 

HEAD - Data File Header 

This is a routine which allows the user to examine and/or 
change the contents of the data file header. For example, changes 
may be made in initial clement (IE), last element (LE) , initial 
line (IL) and/or last line (LL) if necessary, for a specific 
appl icat i.on . 


32 


PMAT - Print Matrix 


This is an ELAS module which prints a matrix of data about 
a selected point from an ELAS data file. As an example, this 
module was used in the accuracy evaluation of the Shenandoah 
classification, and a print-out from this module is shown as 
Figure 6. The size of the matrix is variable and the 5.nput 
for the entire point may be either by element and line number 
or by UTM coordinates. 

Landsat Data DESTRIPE 

The process of destriping should be applied only to data of 
the Goddard format ^diich shows a banding pattern in the raw data 
on one or more bands. This module consists of two overlays, 

DSTB and DSTK. 

COMP - Common Image Display Module 

This module is a general purpose image display module used 
for displaying data and manipulating data for display. It is in 
this module that one may build a color table for a classification 
that has been complete by SRCH (see below) or some other overlay 
or combination thereof. Also, it is possible to slide color in 
this module. The command, SC can be used to highlight successive 
classes in the color table. The highlight color is rolled one 
position forward in the color table of the classification each 
time the carriage is returned. 

In this module, as in most others, typinj in LD on the console, 
will list the directives available in COMD. Some of these are 
briefly reviewed below: 


33 


Figure 6 


ORIGINAL PAGE IS 
OF POOR QUALITY 


30f. 

307 

30 A 

309 

310 

311 
M'd 
313 
3m 
315 

31 A 
317 
31H 

319 

320 
PMAT 


CEDirn 748M75 430537A 

148414651488146714681489147014711472147314741475147614771478 


2 2 2 4 15 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 ^22 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 
2 2 2 2 2 


5 22 6 6 6 
4 4 2 2 15 
2 4 2 2 5 
2 2 2 5 1 
2 2 2 15 1 
2 2 2 2 8 



2 2 2 2 2 

2 2 2 2 2 

2 2 2 2 2 

2 2 2 2 4 

4 4 2 2 4 


6 5 3 4 4 
6 5 4 2 4 
5 6 5 2 4 
2 6 5 2 2 
2 5 5 2 2 
8 2 2 2 2 
8 2 15 6 8 
2 415 2 6 
2 4 4 2 6 
2 2 4 6 1 
2 2 4 6 1 
2 2 4 2 6 
2 3 8 6 6 
<*8866 
4 15 2 2 2 


34 


List directives. 


a . 

b. 


c . 


d. 

e . 


f . 


ft- 


h. 


LD - 

DP - Display picture on image display device. This 
image is displayed through both the function 
and color table currently in memory. 

BT - Build Color Table. This directive permits 
creation of a specific color for a class or 
range of classes. After this directive is given, 
the user is prompted by START, STOP, B, R, G, 
where one will designate the ranges of classes 
to be assigned this color, and the values to be 
assigned to the blue, red and green guns of a 
CRT. This range of values is from 0 to 15. 

SC - Slide color - as described above. 

ST,- Store color table name with a four character 
name , e. g. , CMAT. 

LT,- Load color table which was previously stored with 
a given name . 

PT - Print values assigned for all classes, 0-63 for 

the blue, red and green given values of a specific 
color table. 

DF - Display function. Phis function determines the 
manner in which the color table is displayed 
against the data. Normally, a linear function 
is used with a slope equal to one so that the 
value assigned to a class in the table is also 
equal to the class which is displayed by the 
imaging device. Other options, however, are 
available, sucli as a linear function with a slope 
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i . CK 
i . KF 
k . LF 
1 . KC 

ni. CV 


different fron unity, a five function or an 
arctangent function. 

Compute function. The user may change the 
fvmction through manipulation of the trackball. 
Freer.e function. The function produced is rc- 
taineil in a non-alterable state by this commanii. 
l.oad function. The function \\/hose name is 
stored in the subfile, is loade».i. 

Read coordinate. The pixel under the upper 
left (pK'ulrant of the cursor is read with the 
coordinate given in the header of the data file 
tin KhAS . each data file lias a header record 
which describes the data in the file). 

The cursor moves througli t ite data displayed 

.md stops so that the pixel under the upper left 

qu.ulrant of t lie cursor is at the input coordinates. 


l./\Hl. - h.ibe l_hnogq Cap.ibiJ ij^^ 

riiis iiuHiulo is .1 set of editing and labeling routines for use 
with an imaj’.o displ.iy device. With these routines, it is possible 
to lu'ing exp lan.it ofv a Ijdia -numeric information into .in image. 


Mul t i t em por.i 1 Kogist r^t ion oj" l.ands at Pat a 

I’here .iii' six K1./\.R modules which accomplish scene - 1 o-scene 
rei’istr.it ion i>f 1,. inds.it d.ita. (,)ne set . (OCON. FTIIT .uiil OVl./\) w.is 
written \ov registering data in the Coddarii format . A parallel 
set O’MCO. rMFll aiul I’MOl.) was written for registering dat.i in tlie 
KROS form.it . Onlv t lie first set will be discussed here. 


In 


CKION - Comp ut o M.tpplng Coefficients, Goddard Fommt 

This module Is used to set control points in the dats. A 

control point is a pixel in the data set Identified by scan line 

and element number, whose UTM coordinates can be determined from 
a map. An example might be the intersection of two malor high- 
ways where this intersection is visible in the data set. 

PTU T - Toint Hu nt , doddard Format 

In operation. PTHT is used to aut oma t lea 1 1 y sfloct control 
points for l.andsat scene- 1 o-scene registration. 

OVLA Oy e r One Scen e to Another - dodd^ird Format 

This module will map data from one l.andsat frame into the 
coordituUi' svst em of .'mother. 

SUC H - S earclA 

I’l Kl.AS , SRilH is .m automated procedure for acquiring spec- 
tr.'illv honu\genous trainin.g fields fntm nmltlvariate data by passing 
a 1 X 1 window through the data. The rational behind this program 
was (.1 i scusst'vi c.jrl ier (.Clbula, 1*^801 \^?lth the dat.'i for only the 
.1 X .1 winilow being applical->le to SRCM as found in KIJVS . 

FTF Al. - ^i'o i ny t.' 1 us t y r 

Witli I'rt'Al.. t lu' user may col lect tr.'iining st.it 1st ics by point- 
bv-jHunt clustering within a polygon of the data set . PTCAl. may be 
considered .is .i version of SRCH that allows input of one point at 
.'I t inu> t\'r mergiM'. 

SUP F. - Superv istni Kn cov i n te r 

This is a module in KIAS which will permit the user to develop 
statistics with .i superviscil approach. SUPK computes st.it istics 
im .ill I'lemenis foi’ .ill channels within a defliu'd polvgon (contained 



in subfile HOTP) . This signature may Chen be added to the 
statistics stored in the stat file (STF) . 

STPR - Statistics Print 

STPR is the ELAS module which will permite the user to print 
the statistics produced from SRCH, PTCL or SUPE described above. 

COLR - Color 

This module permits the automatic development of a color 
table based on the statistic*’ stort j in subfile STF. In practice, 
the user is requested to input whi h two channels are to be used 
from subfile STF (with Landsa^ data, usually one channel is from 
the visible region of the sp .trum, while the second is one of 
the two infraf’d channels). In the PLOT option, a two-dimensional 
representation of the data is produced, using the first channel 
as the Y-axis and the second as the X-axis. The class means for 
the inputted channels became the x and y coordinates of this graph 
(this is the two space plot described earlier by Cibula, 1981). 

With this overlay, colors are then assigned on the basis of the 
positions of each of the class means. 

Mi\XL - Maximum l.ikelihood Classification 

With the SRCH, PTCL, and SUPE modules, classes of signatures 
are developed whoso signatures are placed in subfile STF. MAXL will 
use these statistics to classify the data on the basis of maximum 
likelihood. Each pixel examined is assigned to the class where the 
statistics of that class best fit the irradiance values of that pixel. 

Mapping Landsat Data to the Universal Transverse Mercator ( UTM) Grid 
There are four modules which can be used for mapping Landsat 
data to the UTM griil. Two modules (OC.CN and OGEO) arc used with 
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with Landsat data in the Goddard format while the other two mod- 
ules (PMGC and PMGE) are used with data in the EROS Data Center 
format . 

The first module of each format (OGCN and PMGS) employs con- 
trol points and computes Landsat coordinates as functions of 
northings and eastings. The second module of each format (OGEO 
and PMGE) resamples the pixels of the Landsat grid into a disc 
file that represents a UTM grid, with the resampling being done 
on a cell-by-cell basis. 

POLY - Polygon Selection 

In ELAS , POLY is a module which is used in picking polygons 
from an image display device and saving them in HOTP . HOTP is 
a subfile used for storing polygons. The selection of another 
polygon will cause the previous one in HOTP to be overwritten. 

In order to save a polygon initially stored in HOTP for future 
use, it must be stored in subfile PGF. To use a previously 
stored polygon, this polygon needs to be moved from PGF to HOTP. 

PGED - Polygon Edit 

This module is the ELAS polygon edit module and operates 
only on the polygon in HOTP. 

PLYA - Polygon Acreage 

This module determines the frequency of occurrence of a 
particular class or group of classes within a designated polygon 
imbedded in a data set. 

PLYX - Polygon Extraction 

This module is used to extract the contents of a polygon. 
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PGUP - Polygon Update 


This is the ELAS polygon module. This permits defined changes 
to be in the data only within the boundary of a specified poloygon. 
For example, if the user wished to change classes 12, 21 and 28 
to class 75 only within a restricted area, this area wculd be 
enclosed by a polygon and the desired change would be effected 
by application of PGUP. 

Pata Base Analysis and the Programmable Calculator 

In ELAS, a data base is simply a data file containing all 
the components that are required by a user to perform a particular 
application. As an example, a data base may consist of digitized 
elevation data (available from NCIC tapes) which corresponds 
point for point to a completed SRCH Landsat application. Each 
component (the completed SRCH classification and the elevation 
data), exists as a channel of the data file and can therefore, 
be displayed as an image. Continuing with the example above, 
differing elevations could be displayed in color with differing 
shades and hues of color if one had constructed an appropriate 
color table. Furthermore, the data base could contain other 
channels of information such as slope and aspect (derived from 
the NCIC elevation data through application of the T0P6 module) 
if these data are necessary in the calculations for a specific 
application. In addition, it may be necessary to define ir- 
regular polygons that each represent some environmnetal or other 
factor that relates to the desired application. Once the vertices 
of these polygons are determined, the data from the digitizer is 
reformatted and written into the PGF subfile. Then the PGUD module 
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is employed to write this data into one channel of the data 
base. These polygons should lie within the rectangular region 
defined by the data base. In application to the Shenandoah Land- 
sat study, all of these data were required for the final products. 

The normal method of creating the data base file is through 
the use of the GEOREF modules or the TOPO module (discussed below) . 
If UTM boundaries are specified in these modules, the header is 
automatically built. This was the mode employed in this study. 

In most cases, naturally occurring vegetation is determined 
by the complex interaction of soil types, climatic factors, and 
topographic features. A simple example of the application of 
the TBED and PCAL modules in the ELAS subsystem as applied to the 
Shenandoah data would be the breakout from the general grassland 
and pasture class (Class 1), the lowland pastures and grasslands 
(below 2,230' in the valley from Big Meadows and other high ele- 
vation meadow areas (above 2230'). For purposes of our discussion 
here, we will designate these high elevation meadows as Class 53. 

In use, the operator would designate t!ie data base file as input 
(IDI) as well as output (ODF) in the FMGR modulo . With the TBED 
module, the user builds an index table that indexes each value 
to itself except for location 1, which is set to value 53. This 
is stored as Table 1. In the PCAL module, in the SP mode, the 
user would specify the number of channels (register) as 3, des- 
ignate these registers as channels 1, 2, 3 where channel 1 
contains the coded elevation data, channel 2 the Landsat classi- 
fications and channel 3 is specified as being the output channel. 
The table numbers are specified as 0, 0, 1 where an 0 indicates 
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that the corresponding channel will not be indexed through a 

( 

table and with 1 indicating that this channel is indexed through 

table 1. •• 

The rCAL instructions are then: '* j 

LM R1 

- 45 (Value 45 is the code which corresponds to an 
elevation of 2,250 ft) 

IF 4 4 6 " ‘ 

LM R2 

i:x 

LM R3 
EC 

The execut ion of these instructions will cause all cells in 
the Landsat c 1 ass i f i cat ion which arc class 1 and below 2,250 ft. 
to be left unchanged. Those cells above 2,250 ft. that were class 
I are now changed to class 5i. All other classes remain unchanged. 

With this operation complete, one could now designate class 1 as 
lowland pastures and grasslands and class 53 as high elevation 
meadows . 

In a similar manner, other breakouts could be achieved (e.g., 
separate tlie lowland pastures from the high elevation red oak 
forests witli newlv emergent leaves or the coded aspect and slope 
values could be used in the same manner as the elevation data if 
models existed which sliowed the interrelationships between these 
variables . 

TRED - Table Ed itor 

TRED is the module used in ELAS for editing functions and 
indexed tables and has two basic capabilities. The primary 
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function is the editing of indexed tables (EIT) and functions 
(EF) . The non-editing function of TBED if format conversion. 

One command (ITTF) will cause an index table to be switched to 
a function and the other command FTIT, causes the inverse. 

PCAL - Programmable Calculator 

This module allows the user to program the desired algorithm 
at run time. In use, an instruction set is input to define the 
function or functions to be processed. Each pixel is processed 
according to these instructions (e.g., "If a given pixel is class 
1 and the corresponding pixel in the elevation channel shows that 
this elevation is above 2,250 feet, change this pixel to class 53") 
and an output is made for each cell. 

NCIC - Topographic Data Manipulations 

Three modules comprise this aspect of ELAS which will permit 
the construction of topographic data bases. The first module con- 
tains six overlays, TOPO , TOPI, TOP2, T0P3 , TOP4 , and TOP5 . The 
output from the T0P5 overlay is a channel of data which contains 
pixels whose centerpoint corresponds to specific coordinates of a 
UTM grid. These pixels may be of variable size, depending upon 
the requirements of the data base (in this study, pixel size was 
50M X 50M) . Each pixel is assigned a level value to it which 
corresponds to the mean elevation range over the range of pixel 
size. T0P5 serves as input to the T0P6 module which computes 
slope, aspect and slope length. A third module T6CH computes 
the average north-south slope and the average east-west slope 
in addition to the four variables computed in T0P6 . 
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Digitized elevation data are obtained from the National 
Cartographic Information Center (NCIC) and these tapes serve as 
input to the TOPO series of modules. These modules reformat 
the NCIC data from tape into the ELAS data file format , rotate 
the data 90° to a north-south orientation, compute mapping 
coefficients that relate transposed plate coordinates to UTM 
coordinates (northings and eastings) and resample the 16-bit 
data to a UTM grid file. This file is then used to compute 
slope, aspect, slope length and average slope. Figure 7 shows 
the module flow for TOPI - TOP6 . 

TOPO - Topographic Data to UTM Grid 

This module is the controller for overlays TOPI - TOPS. 

The input is the digitized NCIC terrain tape , and it is in this 
module that the final output cell size is chosen. 

TOPI - Reformat 

This overlay reads the NCIC tape and writes the results 
into a disc file INFl preserving row and column coordinates 
on the tape. 

TC?2 - Rotate 

In this subroutine, the data stored in INFl is rotated 
to a north-south orientation and is stored in INF2 . 

T0P3 - Computer Mapping Coefficients 
This is the subroutine which determines the mapping 
coefficients that relate the transposed plate coordinates 
of INF2 to UTM grid coordinates. 

T0P4 and T0P5 - Resample to UTM Grid 

These overlays resample the data in INF2 to the UTM 

GRID. The output from T0P5 can be visually displayed 
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with an appropriate color table on an image display 
system as a channel of data in a data base. This out- 
put is also used as the input for T0P6. 

T0P6 - Topographic Data to Slope, Aspect and Slope Length 

This routine uses the input from TOPS to calculate slope, 
aspect and slope length. Each of these outputs can also be 
displayed on an image display device as a channel of data in a 
data base. The standard output from this routine is a 4-channel 
file of data for display, classification enhancement (use of PCAL 
option) or for input to other programs. To process a file of 
approximately 4 million cells through T0P6 to produce a 4-channel 
digitized output of elevation, slope, aspect and slope length, 
requires approximat le y 1 hour and 25 minutes of wall Lime on the 
Perkin-Elmer 8-32. 


T6CH - Topographic Data to Elevation, Slope, Aspect, Slope Length , 
Overage N-S slope a. id Overage E-W slop e 

This routine is a modified version of T0P6 and uses the output 
from T0P5 to produce an output file of data containing the following 
six channels: 

Channel 1 - Elevation 

2 - Slope 

3 - Aspect 

4 - Slope length 

5 - Average north-south slope 

6 - Average east -west slope 
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APPENDIX II 


CRiGMMAL PAGE IS 
OF POOR QUALITY 


As discussed on page 19, the first PCAL operation involved 
a situation where the elevation data was juxtaposit’ oned over 
the summer SRCH classification. This first calculation involved 
several elevational breaks where differing land cover types at 
differing elevations were included within the same spectal class. 
These spectral classes represented pastures at low elevations 
while above 2500 ft. either represented high elevation meadows 
or red oak communities. The first calculation was performed 
on classes which were definately known to be in need of an eleva- 
tional separation as well as several which were chosen where it 
was desired to examine the resultant separation in mr e d* 'ail. 
These separations are shown in Table VI . 


Table VI 


Original 

Spectral Class 

8 

11 

15 

16 

17 

18 

19 

20 

22 

24 

25 

26 

27 

30 

31 

New Class 
Below 2500' 

8 

11 

15 

16 

17 

18 

19 

20 

22 

24 

25 

26 

27 

30 

31 

New Class 
Above 2500' 

36 


39 

40 

41 

1 

43 

44 

45 

46 

47 

48 

49 

50 

51 


With the previous operation complete, it was then decid<^d to 
address two discrete areas within the park which had suffered defo- 
liation probably is a result of insect damage. Within polygons 
which defined these areas, completion of the PCAL operation yielded 
class 52 which represented this defoliation which earlier was 
represented by a portion of spectral classes 23, 28 & 29. This 
was accomplished by the application of overlay POLY and PGUD. 
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At this point, an examination of the winter data was \mder- 
taken. Very few spectral classes in the winter frame appeared to 
have a high ratio of IR to visible reflectance. Those spectral 
classes which did, and were related to areas within the park were: 

1, 6, 21, 24, & 26, Field studies designed specifically to ad- 
ress these classes were completed and the results of this analysis 
are given in Table VII. 

Table VII 

FIELD ANALYSIS OF SPECTRAL CLASSES FROM FRAME 2401-15102 

WHICH HAVE A HIGH IR/VIS REFLECTION RATIO 

Class Description 

1 Hardwoods with evergreen understory, primarily red 

oak with mount ian laurel ( Kalmia latifolia ) 

6 Denuded pastures, moderate to light inert materials 

21 Dense pine ; some areas mixed pine/hardwood but high 

percentage of pine in all cases 

24 Hardwood canopy (primarily red oak) with an understory 

of mountain laurel 

26 Mixed pine/hardwoods 

Field data from both the summer and winter data sets were 
analyzed to provide input for various PCAL operations where several 
different logical strategies were used. For these operations, the 
winter data was superimposed over the summer data set so that there 
was a one to one correspondence for cells in each data set that 
shared the same UTM coordinates. In one case where elements of a 
certain group of spectral classes from the summer data (Set A) 
superimposed with elements of a particular soectral class in the 
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winter data, this combination was transformed into a previously 
non-existant class. When this same winter spectral class was 
coincident with yet another set of classes (Set B) of the summer 
data, this second combination was changed into yet another non- 
existant class. There was a third set of spectral classes which 
could be designated in the summer data iSet C) , where when other 
winter spectral class coincided with any of these, no change was 
effected. These operations are detailed in Table VIII (Set C is 
implicit in this table) . 

With this operation completed, a partition was m.ade with 
spectral class one. Below an elevation of 2250', this class 
remained one, above this elevation, this class became 53. This 
separation helped to better define Big Meadows and other high 
elevation meadow areas. 

As a result of these operations, there existed a total of 
63 computer classes. This is the maximum number of classes which 
could be accommodated on our system. At this point, application 
of ground truth data permitted the 63 total classes to be grouped 
in such a manner that a new classification of 32 grouped classes 
remained. All subsequent operations were performed on this 
grouping where the grouped classes 1 through 32 were considered 
as mathematically operable classes as was the case with the 
spectral classes discussed above. These new groupings are shown 
, in Table IX. 

With this new ■'gioc.pir.g ■ an' inp‘ut , grouped class 32 
(representative of spectral classes 28 & 29) was left as class 
32 below 2000'. Above this elevation, it was included in the 
aggregate of grouping 19. This high elevation break continued 
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TABLE VIII 


ORiUlNAL PAGE IS 
OF POOR QUALITY 


MANIPULA.TION FOR VARIOUS PCAL OPERATIONS 


PCAL 

Operation Description 

la If winter spectral class 1 associates with summer classes 

2.3.4.5.6.7.10.12.13.32.39.43.46 to 49 or 53, these com- 
binations become class 54 

lb If winter class 1 associates with any other class except 

class 52, then these combinations become class 55 

2a If winter class 21 associates with 2,3,4,5,6,7,10,12,13, 

32.39.43.46 to 49 or 53, then these combinations become 
class 58 

2b If winter class 21 associates with any other class except 

class 52 or 53, then the combinations become class 59 

3a If winter class 24 associates with 2 to 7 , 10, 12, 13, 32, 

39, 43, 46 to 49 or 53, then these combinations become 
class 60 

3b If this class associates with any other except for classes 

52, 55, and 59, then these combinations become class 61 

4a If winter class 26 associates with 2 to 7, 10, 12, 13, 32, 

39, 43, 46 to 49 or 53 then these combinations become 
class 62 

4b If this class associates with any other except for 52, 55, 

59 and 61, then these combinations become class 63 

5 A transformation of spectral classes 7 and 12 as shown 

below: 


Input classes 

7 

12 





Final Class 

f 7 

12 


Assignments 

J 


1550' 


1 38 

38 

T50t5 




6a If winter class 6 associates with 1, 8, 10, 11, 15-17. 

19, 23, 28-31 then tricse combinations oecome class 56 

6h If winter class 6 associates with 9, 14, 34, or 35, then 

these combinations becon.e class 57 



TABLE IX 


Grouping From 63 Classes to 32 Groups 
(Intermediate Classification) 


Group 

Classes Included 

Description 

1 

1 

Grass & forbs below level 45 

O 

2.3.5.43.46.47.48. 

49 

Red Oak forests 

3 

4.13 

Yellow poplar; oak/hickory 

4 

6.10.21 

Chestnut oak/red oak 

5 

7 

Red oak/chestnut oak 

6 

8 

Old field succession 

7 

0.14 

Appears to be red oak/chestnut oak 
mix; mav include some apple orchard 

8 

11.30.31 

Dark inert 

q 

12 

Red oak/chestnut oak 

10 

15 

Very small amount - appears to he 
red oak 

LI 

16.35 

Moderate inert 

12 

17.19 

Moderate, bright inert 

13 

18 

Pine 

14 

20.22,24.25.26.27 

Red oak/various mixtures 

15 

23 

Bright inert 

16 

32.30 

Red oak forests 

17 

33 

Vigorous forest regeneration; some 
yellow poplar 

18 

34 

Very dark inert 

19 

16.51.53 

Upland meadows - Big Meadows 

20 

37.38.40.41.42.44. 
4 5.50 

Red oak 

21 

52 

Defol iat ion 
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Group 

Classes Included 

Description 

22 

54 

Hardwood (R.O.) with evergreen under- 
story 

23 

55 

R.O. with evergreen mix or understory 
possible grass understory 

24 

56 

Inert with slight vegetation present to 
pastures: definite lowland class only 

25 

57 

Pasture 

26 

58 

Dense pine 

27 

59 

Mixed pine /hardwood 

28 

60 

Hardwood evergreen understory; KALA 
or FIST 

29 

61 

Pine or pine/hardwood mix 

30 

62 

Pine to pine /hardwood mix 

31 

63 

P ine 

32 

28 & 29 

Lowland inert or meadows (below 2000’') 


1 

I 
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to develop the concept of high elevation meadows. The extent 
of these meadowlands was not yet satisfactory, but was solved 
by the use again, of the winter data set. In this data set, 
spectral classes, 4, 6, 8, 10, and 20 related to the upland 
meadows above 2000'. Consequently, the following logic was 
used: If classes 4, 6, 8, 10, or 20 above 2000' elevation are 

over any class in the composite summer data except classes 8, 11, 
12, 15, 18, or 19, make these occurrences class 33. This oper- 
ation served to give a separation between those areas of upland 
meadows which were essentially mixed grasses and forbs (19) from 
those areas which were mixed grasses and forbs as well as recently 
emergent perennial vegetation. 

With these operations complete, analysis of field data and 
this grouping using the option of ELAS , a new grouping of 22 
combined classes was performed. During the course of this 
final analysis, it was recognized that spectral class 8 which 
was in grouping 5 was significantly different from spectral 
classes 17 6t 19, so that a separation was needed. Therefore, 
the original data was retrieved and spectral class 8 was placed 
in grouping 22 (old field succession) in the final classification. 

The final groupings of this classification are given in 
Table III and color coded classification of these final land 
cover types is given with the legend in Figure 5A and 5B. 
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Appendix lii 

Field Sites Derived From SRCH Classifications 

Summer Sites 


UTM 


Sl-03 

702100 

Sl-50 

724200 

Sl-105 

723250 

Sl-156 

738450 

S2-04 

•^01200 

S2-05 

700750 

S2-06 

706050 

S2-51 

718550 

52-52 

716450 

52-53 

722900 

52-lOt. 

736950 

52-107 

735200 

52-ioa 

741550 

52-156 

749000 

52-157 

747450 

53-07 

705150 

53-08 

704700 

53-09 

716500 

53-54 

711300 

53-55 

71530U 

53-58 

710750 

53-109 

733000 

53-110 

753600 

53-158 

747150 

54-10 

696400 

54-11 

698650 

54-12 

702750 

54-57 

702650 

54-58 

715350 

54-59 

7167L0 

54-111 

7325‘=>0 

54-112 

732600 

54-113 

732450 

54-114 

732600 

54-159 

746450 

54-160 

750350 

55-13 

699900 

55-14 

701450 

55-15 

716750 

55-60 

721000 

55-61 

723050 

55-62 

722250 

55-115 

733400 

55-116 

740350 

55-117 

726400 

55-161 

7449PC 

55-163 

757350 

55-164 

759350 

56-16 

697700 

56-17 

700700 

56-63 

715300 

56-64 

708900 

56-118 

732650 

Sb-1 19 

751100 


42iJ95bO 20 15 

42‘+5950 

4285750 

42^9600 

42197C0 

42iJ5550 

428505U 

4247d50 

4246300 

4Pf»055f! 

4283450 

42^4400 

42*53050 

4303200 

4302700 

42<i6700 

4233000 

4c35100 

42*+b250 

42*+8a50 

4252000 

4282000 

4282100 

4303000 

42*4150 

4227700 

4235600 

4244150 

4251150 

4252650 

4282900 

428225U 

4289800 

4290450 

4292700 

4306o50 

42l9bCn 

42*^4900 

4234800 

4248700 

4242800 

426C450 

4289950 

429400C 

4279600 

4294700 

4306150 

43U44C0 

42-^5200 

4234650 

4246250 

4241600 

4284350 

4291300 


ORIGINAL PAGL 
OF POOR QUALITY 
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Sb-165 

Sfa - 1 66 

S7-18 

S7-19 

57-65 

57-66 

57-120 

57-121 

57-167 

57- 166 

58- 19 

58-67 

58-b8 

58-122 

58-123 

58- 169 

59- 20 

59-69 

59-124 

59-169 

59-170 

510-21 

510-22 

510-70 

510-71 
510-125 

510- 171 

511- 23 

511-72 

511-12b 

511- 127 

512- 24 

512-25 

512-73 

512-74 
512-128 
512-129 

512- 173 

513- 2b 

513-27 

513-75 

513-130 
513-131 
513-132 
513-174 
513-175 

513- 176 

514- 28 

514-76 

514-133 

514- 177 

515- 29 

515-77 

515-78 

515-134 
51b-30 


751350 

750650 

695700 

702200 

711950 

717900 

742350 

735750 

742350 

740750 

706600 

714550 

711150 

722200 

722500 

747050 

697950 

727500 

737400 

750750 

743250 

697900 

696650 

717250 

717250 

729350 

753350 

707950 

701950 

732500 

744100 

703200 

703150 

718450 

714350 

736150 

735900 

735150 

700900 

700150 

723750 

730900 

739300 

742600 

757350 

758200 

743000 

700900 

720000 

738600 

748100 

695000 

712100 

723200 

734500 

694700 


4304500 

4306950 

4225400 

4212150 

4280200 

42b0050 

42^4850 

4276300 

42^4850 

42^4500 

4228350 

42^0800 

42b7150 

4265150 

42^0250 

4314100 

4219450 

42^7550 

4292900 

4310850 

4295C00 

4226950 

42-30000 

4256650 

4258200 

4295450 

430710U 

4219650 

4249300 

4278450 

4283750 

4234750 

4220250 

4253800 

4246250 

4294400 

4277250 

4293400 

4231100 

4228050 

4256700 

4289250 

4293000 

42900CO 

4304600 

4313550 

4305500 

4232100 

4244200 

4284250 

4298600 

4225050 

4262250 

4263950 

4292100 

4223500 
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S16-79 

516- 135 

517- 31 
S17-80 
Sl7-13fa 

817-173 

817- 179 
816-32 
816-61 
816-137 

818- 180 

818-181 

819-33 

819-82 

819-63 

819-138 
819-182 

819- 183 

820- 3«4 

820-84 

820-139 

820- 184 

821- 35 

821-85 
S21-140 

821- 185 

821-186 

822- 36 

822-86 

822-141 

822- 187 
822-188 

823- 01 

823-02 

823-87 

823- 169 

824- 37 

824-68 

824-142 

824-190 

824- 191 

825- 89 

825-90 

825-143 

825-192 

825- 193 

826- 38 

826-91 

826-92 

826-144 

826- 194 

827- 39 

827-40 

827-93 

827-145 

827-195 


721000 

730750 

707650 

708550 

732300 

749100 

746300 

704250 

706100 

742600 

757100 

758500 

711750 

707350 

704200 

727000 

746050 

736100 

699400 

703550 

743000 

743100 

705250 

718950 

732100 

740750 

754950 

698100 

714350 

733100 

743250 

758400 

703900 

710700 

704300 

736700 

694750 

717100 

734850 

742100 

736900 

713900 

722900 

727250 

752150 

742450 

691300 

717550 

724850 

728500 

743650 

697850 

705350 

717900 

727250 

742450 


4258250 

4274350 

4215350 

4255500 

4278250 

4311050 

4306900 

4214400 

4252050 

4279300 

4301650 

4302400 

4227100 

4258200 

4250300 

4270250 

4311450 

4302600 

4232300 

4253950 

4292500 

4295300 

4230150 

4244100 

4288650 

4296200 

4306600 

4232050 

4245600 

4279550 

4295550 

4294800 

4218700 

4229200 

4251400 

4304300 

4232950 

4252100 

4292150 

4295450 

4292700 

4245550 

4263050 

4274500 

4296450 

4310400 

4224350 

4249400 

4261600 

4273250 

4296100 

4223250 

4216600 

4249050 

4280000 

4295100 
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S26-41 

S26-94 

S28-146 

528- 19b 

529- 42 

S29-43 

S29-95 
S29-147 

529- 197 
S30«96 

530- 97 

S30-148 

530- 198 

531- 44 

S31-98 

S31-99 

S31-149 

531- 199 

532- 45 

S32-4b 

S32-100 

S32-101 
S32-150 
S32-151 
S32-200 

532- 202 

533- 47 

S33-102 

S33-152 

533- 2U1 

534- 48 

S34-103 

S34-lh3 

534- 203 
S3b-49 
S3b-104 

535- lb4 

S35-204 
»XPR 


704150 

721750 

728500 

752500 

698900 

704250 

706500 

741850 

744250 

718450 

725000 

729500 

740900 

703000 

702750 

716400 

740000 

744600 

692950 

696600 

720750 

721250 

729250 

734450 

745500 

738250 

696750 

723150 

722750 

735150 

7U6U50 

712400 

739550 

735700 

702000 

716900 

733950 

744360 


4216000 

4247100 

4294700 

4298050 

4214850 

4217850 

4251900 

4283250 

4311700 

4253100 

4264400 

4274350 

4299300 

4230300 

4248750 

4262100 

4276700 

4313400 

4227350 

4230350 

4248950 

4256650 

4271650 

4275850 

4296950 

4293200 

4224150 

4254400 

427545C 

4297350 

4232350 

4261200 

4282900 

4312400 

4233900 

4255300 

4292350 

4300400 
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Fall Sites 
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fall 

SITFS 



F05-01 

4214] 00 

685200 

WAYNESBORO WEST 

F05-18 

4244250 

697350 

rCGAHEYSVILLE 

F05-30 

4238400 

722000 

standakdsville 

F05-32 

4270100 

7193J0 

niG MEADOWS 

F08-17 

4219b00 

701600 

crozet 

F08-39 

4285100 

713600 

STANLEY 

F08-82 

4270750 

723800 

BIG MEADOWS 

F08-63 

4278250 

733150 

OLD RAG MTN. 

F12-07 

4218000 

695100 

WAYNESBORO EaST 

F12-22 

4245801.' 

704550 

mcgaheysville 

Fl?-48 

4271200 

729400 

OLD RAG MTN, 

F12-49 

4270550 

729000 

old rag MTN. 

F13-08 

4228950 

694950 

CPIMORA 

F13-23 

4239950 

698200 

mcgaheysville 

F13-36 

4284000 

714250 

STANLEY 

F13-50 

4277750 

729450 

OLD rag MTN, 

F14-02 

4223400 

695250 

CRIMORA 

F14-19 

4232300 

698950 

BROWN'S COVE 

F14-31 

4237050 

715100 

SwiFT RUN GAP 

F14.35 

4266750 

722800 

DIG MEADOWS 

F17-03 

4221350 

690300 

WAYNESBOROiEaST 

F17-20 

4232650 

697450 

BROWN'S COVE 

F17-34 

4271300 

729000 

OlO rag MTN. 

F18-09 

4228500 

690800 

CRIMORA 

Fla-24 

4244500 

701200 

mcgaheysville 

Fl8-bl 

4272300 

728300 

DIG MEADOWS 

F19-06 

4223150 

697950 

BROWN'S COVE 

F19-45 

4279950 

727250 

LJRAY 

F19-48 

4278350 

728200 

luray 

F19-47 

4276550 

729100 

thorntom gap 
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F20-12 

4230350 

695450 

CRIMORA 

F20-26 

4240300 

703100 

mcgaheysvillf 

F21-13 

4227450 

694450 

CRI80RA 

F21-53 

4290450 

731100 

thornotn gap 

F22-27 

4252700 

708700 

ELKTON EAST 

F22-5H 

4275000 

721550 

biG meadows 

F23-14 

4220450 

698750 

CROZET 

F23-55 

4287200 

728100 

luray 

F25-10 

4223600 

694450 

CRIMORA 

F25-11 

4227600 

695600 

CRIMORA 

F25-25 

4239950 

704800 

mcgaheysville 

F25-36 

426240C 

719050 

FLETCHER 

F25-52 

4290850 

730850 

THORNTON GAP 

F26-04 

4224000 

704600 

BROWN'S COVE 

F26-05 

4227650 

703100 

BROWN'S COVE 

F2f,-3ri 

4267850 

725650 

biG meadows 

F26-41 

4290150 

726^00 

luray 

F26-4? 

4288600 

731200 

THORNTON GAP 

F2fa-43 

4276800 

730000 

OLD RAG MTN, 

F26-44 

4278350 

732550 

OLD RAG MTN. 

F27-15 

4229000 

697450 

BROWN’S COVE 

F27-5t 

4276450 

731100 

old rag MTN. 

F27-57 

4284400 

7^4150 

THORNTON GAP 

F27-56 

4284500 

7 53700 

THORNTON GAP 

F28-lf> 

4219900 

635750 

WAYNESBORO EAST 

F28-28 

4245950 

704200 

mcgaheysville 

F28-2S# 

4248400 

708550 

SWIFT RUN GAP 

F28-59 

4274650 

730650 

OLD RAG MTN. 

F26-b0 

4290500 

732250 

THORNTON GAP 

F28-61 

4290200 

730400 

THORNTON GAP 
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WOI-74 

4244000 

714150 

SWIFT KU^i G8H 

WOl-63 

4294300 

740700 

CHESTER GAM 

U02»69 

4229fa!j0 

692750 

CRIMORA 

W02-75 

4260200 

712650 

ELKTOM EAST 

W02-84 

4289150 

734650 

THORNTON GAP 

W04.76 

4248100 

702550 

MCGAHEYSVlLLt 

H04-85 

4279050 

724000 

LIJRAY 

W05-77 

4263850 

723150 

fletcher 

W05-8b 

4280200 

732250 

THORNTON GAP 

W06-76 

4253350 

705050 

ELKTON WEST 

WOf-87 

4280000 

744650 

WASHINGTON 

W21.64 

422l8u0 

697300 

CROZCT 

W2l-b6 

4224700 

700150 

RROWN'S COVE 

W21-70 

4248700 

709850 

SWIFT RUN GA^> 

W21-71 

4249500 

709300 

SWIFT RUN GAP 

W21-79 

4281200 

730850 

THORNTON GAP 

W24-65 

4222800 

695550 

WAYNESHORO EAST 

U24-b6 

4214950 

b96650 

WAYNESBORO FAST 

U’24-72 

42bl lUO 

7159Q0 

ELKTON EAST 

U24.80 

4276200 

730350 

OLD RAG MTN, 

W26-b7 

4224350 

702350 

BROWN’S COVE 

W2fe-73 

4246900 

704050 

mcgaheysville 

W2b-bl 

4273350 

728650 

BIG meadows 

W26-8? 

4278250 

734600 

OLD RAG BTN. 
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Spring Sites 


spriug 

SITCS 



SPl-89 

4229800 

701200 

FROWN'S COVE 

SPl-99 

4265400 

727050 

HIG MEADOWS 

SPl-107 

4291000 

737350 

THORNTON GAP 

SP2-94 

4243250 

700800 

MCGAHEYSVILLC 

SPl2-9f> 

4237700 

695200 

GROTTES 

SP12-103 

4269600 

723600 

BIG MEADOWS 

SPl&-96> 

4234850 

699650 

BROWN'S COVE 

SP16-1 04 

4265900 

718700 

BIG MEADOWS 

SP26-97 

4228650 

696400 

BROWN'S COVE 

SP26.-105 

4261150 

718650 

FLETCHER 

SP27-109 

4290350 

730950 

THORNTON GAP 

SP28-92 

4235000 

704200 

BROWN'S COVE 

SP31-90 

4244450 

706100 

MCGAHEYSVILLL 

SP31-10P 

4264150 

716650 

STANLEY 

SF31-108 

4292050 

738950 

THORNTON GAP 

SP33-93 

4241000 

703500 

MCGAHEYSVILLC 

SP33-101 

426300 ' 

720800 

FLETCHLR 

SP3‘j-88 

4237'ir 0 

7 J 7500 

MCGAHEYSVILLC 

SP35-98 

425680V; 

723550 

FLETCHER 

SP35-10b 

4298500 

735250 

bentonville 

SP3fa-91 

4237500 

706300 

mcgaheysville 
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:! I 

• $ 

FL-?0 

430072U 

742940 

• - I 
1 

FL-31 

R290750 

738010 


FL-32 

42693BU 

736980 

1 

FL-33 

4271760 

729460 

1 

FL-7«4 

4266430 

727760 

t 

FL-35 

4263960 

723060 

• ♦ 1 

FL-36 

4272660 

7281^0 


FL-37 

4272670 

728600 


FL-38 

4265660 

723290 



I 
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Appendix IV 

Random Sites for Accuracy Verification 


1 

696725 

9217376 

2 

722675 

9265626 

3 

75207b 

9312376 

4 

702125 

9295276 

5 

712575 

9252376 

6 

739125 

9285876 

7 

729025 

9251776 

8 

736525 

9289026 

9 

717725 

9251676 

10 

639925 

9217326 

11 

739725 

9287976 

12 

725275 

9259976 

13 

718175 

4292376 

14 

729325 

9280926 

15 

740325 

9298326 

16 

722525 

9268726 

17 

733825 

9280976 

18 

796825 

9292976 

19 

725675 

9273926 

20 

7UP925 

9291926 

21 

713925 

9236126 

22 

706975 

9250776 

23 

697775 

9218926 

24 

733925 

9259926 

25 

738025 

9296326 

26 

713175 

9251626 

27 

721175 

9261776 

28 

731625 

9257226 

29 

713775 

9250826 

30 

716325 

4267876 

31 

708225 

9237626 

32 

710975 

9237676 

33 

710325 

4290726 

34 

7 0 36 7 5 

9225976 

35 

795125 

9303776 

36 

739925 

9282076 

37 

712275 

9258576 

38 

738075 

9286676 

39 

726625 

9259026 

40 

698975 

9290626 

41 

702775 

9232926 

42 

730675 

9265926 

43 

696875 

9237526 

44 

729775 

9259526 

45 

709125 

9292276 

46 

737325 

4,278976 

47 

7 1 9 0 7 5 

92b62'^6 

48 

7216'-5 

4252976 

49 

698225 

9232876 

50 

719525 

9260976 

51 

701925 

9225926 

52 

699025 

4226626 

53 

716625 

4250070 

54 

725925 

925637b 

55 

728675 

4281176 


i 

• i 


56 

696425 

4223576 



57 

710875 

4263176 


' • 

58 

713325 

4236926 
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• » 

59 

60 

734575 

714375 

4268226 

4250076 
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704575 

4242226 



62 

72^075 

4266426 


• ♦ 

6 3 

715125 

4269426 



62 

736375 

4265926 



63 

711675 

4247726 


• « 

6A 

729275 

4292026 



65 

730825 

4261076 


-- 

66 

691425 

4221176 



67 

695075 

4228176 


• • 

68 

735825 

4251876 



69 

7 06675 

4245326 



70 

732225 

4256476 


•• » 

71 

739125 

4299876 



72 

734275 

4282776 


•• 

73 

717825 

4260176 


e • 

74 

692075 

4220426 



75 

711375 

u2u3926 


*. - 

76 

731725 

4277326 



77 

701675 

4242976 


- ' 

78 

739875 

4284376 



79 

730975 

4269576 



80 

712975 

4266326 



81 

711825 

4244626 



82 

733575 

4251826 



83 

732975 

4264226 



84 

746175 

4299976 



85 

704575 

4230626 



86 

735525 

4269626 



87 

734725 

4285126 



88 

696725 

4240576 



89 

740925 

4297576 



90 

715475 

4240026 



91 

728975 

4286576 



92 

702025 

423677b 



93 

734175 

4251076 



94 

739425 

U293676 


• 

95 

710925 

4251576 



96 

716175 

4256376 



97 

693575 

422432b 



98 

7190 75 

4258626 



99 

706225 

4243026 


- * 

100 

737025 

4273526 



101 

745475 

4286026 



102 

696425 

4235176 



103 

711075 

4248476 



104 

731725 

4265726 



105 

754175 

4315476 



106 

712275 

4235376 



107 

733075 

4284326 



108 

720575 

4250926 



109 

718625 

4256326 



no 

718275 

4262476 
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Ill 

736375 

4297526 

112 

693225 

4244476 

113 

749075 

4293026 

114 

751325 

4316226 

115 

748475 

4305376 

116 

727175 

4277276 

117 

717425 

4246226 

118 

704725 

4239126 

119 

751025 

4310326 

120 

732525 

4250276 

121 

706075 

4234526 

122 

726625 

4266476 

123 

609325 

4218076 

124 

726775 

4251776 

125 

741975 

U299126 

126 

71<4525 

4258576 

127 

713025 

4243076 

128 

721775 

4260976 

129 

737775 

4201276 

130 

734925 

4270426 

131 

709675 

4264726 

132 

739275 

4205176 

133 

743925 

4293726 

134 

714825 

4252376 

135 

732475 

4273476 

136 

742675 

4292176 

137 

721325 

4247076 

139 

744825 

U298376 

140 

752225 

4309276 

141 

731425 

4203526 

142 

727775 

4200126 

143 

725275 

4247876 

144 

729625 

4251026 

145 

707875 

4243826 

146 

692325 

U225826 

147 

727225 

4265676 

148 

730675 

4249476 

149 

697175 

4242926 

150 

720425 

4242426 

151 

731425 

4260326 

152 

693675 

4232826 

153 

712375 

4257776 

154 

697475 

4236726 

155 

698225 

4221276 

156 

708925 

4245376 

157 

739275 

4273576 

158 

729175 

4260276 

159 

714675 

4255476 

160 

706675 

4233726 

161 

720675 

4271026 

162 

707725 

4246926 

163 

724025 

4261026 

164 

715425 

4263226 

165 

717825 

4271776 

166 

729775 

4282726 

:67 

743225 

4265976 

168 

750725 

4305426 


169 

690975 

9218876 

170 

716775 

9270226 

171 

703375 

9220576 

172 

729175 

4c?9867b 

173 

732975 

9296676 

174 

693875 

9218126 

175 

708775 

4298976 

176 

692175 

9228926 

177 

740775 

42891 26 

178 

720875 

9267926 

179 

705175 

9291976 

180 

732325 

U276576 

181 

791525 

9296776 

182 

709575 

429962b 

183 

733375 

9278126 

184 

719325 

4269076 

185 

791225 

4279776 

186 

719025 

9293176 

187 

794275 

9287526 

188 

709225 

9250776 

189 

702025 

9225176 

190 

725725 

9250176 

191 

727525 

9297876 

192 

697325 

9259826 

193 

720925 

9259026 

194 

716'^75 

9293926 

195 

751775 

9306976 

196 

797875 

4306126 

197 

730375 

9270376 

198 

718875 

9261726 

199 

709925 

SXFR 

9236076 
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